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AB8TBACT 


The  Naval  Research  Laboratory's  Far  Ultraviolet  Cameras 
Experiment  is  part  of  an  ongoing  project  to  develop  a  means 
for  real  time  space  based  monitoring  of  the  upper  atmosphere. 
The  purpose  of  this  paper  is  to  develop  a  method  by  which 
useable  data  can  be  extracted  from  the  images  obtained  by  the 
cameras.  The  required  output  is  an  intensity  profile  of  the 
light  in  the  band  of  the  camera  as  a  function  of  the  observer 
zenith  angle.  From  this  data  future  analysis  can  be  done  to 
determine  the  spatial  distribution  and  density  of  the 
constituent  species  of  the  upper  atmosphere.  The  long  term 
application  of  this  project  is  to  develop  a  method  for  using 
space  based  imaging  of  the  ionosphere  to  provide  real  time 
radio  frequency  propagation  analysis. 
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Z.  ZMTBOOUCTIOM 


A.  BACXOBOUMD 

The  purpose  of  this  thesis  is  to  develop  a  aeans  for 
analyzing  images  obtained  by  two  far  ultraviolet  cameras 
developed  by  the  Naval  Research  Laboratory.  The  Far 
Ultraviolet  Cameras  Experiment  was  part  of  the  Air  Force 
Program  675  (AFP-675)  payload  flown  aboard  space  shuttle 
mission  STS-39  between  28  April  and  6  May  1991.  AFP-675 
consisted  of  the  following  instruments: 

1)  CIRRIS  >  Cryogenic  Infrared  Radiance  Instrumentation 
for  Shuttle.  This  instrument  measured  spatial,  spectral, 
and  temporal  variations  in  the  Earth's  limb  and  aurora 
between  2.5  and  25 fm. 

2)  FAR  UV  -  Far  Ultraviolet  Cameras.  Designed  to  image 
near  Earth  and  celestial  ultraviolet  sources,  including 
airglow,  aurora,  stars,  nebulae,  shuttle  glow,  shuttle 
thruster  firings,  and  shuttle  H2O  discharges.  This  paper 
will  address  the  airglow  portion  of  this  experiment. 
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3)  URA  -  Unifonsly  Redundant  Array.  Deeonstrated  coded 
aperture  iaaging  of  x-ray  eeiesione  froe  celestial 
targets . 

4)  HUP  -  Horizon  Ultraviolet  Progras.  Measures  spatial 
and  spectral  ultraviolet  character is it ics  of  the  Earth's 
horizon. 

5)  QINMS  -  Quadrupole  Ion  and  Neutral  Mass  Spectrometer. 
Monitors  contamination  in  the  shuttle  environment  and 
provides  input  to  CIRRIS.  (STS-39  PID) 

B.  FAR  ULTRAVIOLET  CAMERAS  BZPERIMBMT 

The  objectives  of  the  Far  Ultraviolet  Cameras  Experiment 
included;  far  ultraviolet  (1050-2000A)  measurements  of  the 
Earth's  upper  atmosphere,  celestial  bodies  including  diffuse 
background  sources,  and  the  shuttle  environment.  The  shuttle 
enviroment  study  included  shuttle  glow,  and  thruster  firings. 
Upper  atmosphere  observations  include  night  horizon  alrglow 
and  the  auroras.  The  purpose  of  these  observations  was  to 
demonstrate  a  means  for  remote  sensing  of  atomic  oxygen  and  NO 
density  versus  altitude  in  the  upper  atmosphere.  The  night 
airglow  observations  were  specificaly  designed  to  demonstrate 
the  ability  to  conduct  remote  sensing  of  the  night  ionosphere. 
The  results  of  this  experiment  will  be  used  to  support  future 
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space  based  real  tise  resote  sensing  of  the  night  ionosphere. 
This  will  be  attespted  by  the  Global  Imaging  Monitor  of  the 
Ionosphere  (6IMI)  project  scheduled  to  be  deployed  aboard  the 
ARGOS  P-91  satellite.  GIMI  will  use  similar  cameras  but  will 
utilize  CCD  detectors  instead  of  electrographic  film  to 
provide  a  real  time  digitized  RF  output.  This  output  will  be 
processed  at  a  ground  station  to  determine  an  altitude 
profile  of  the  atmosphere.  The  work  done  in  this  thesis  will 
provide  the  basis  for  deriving  altitude  and  intensity  from  the 
photographs  and  will  provide  the  background  material  necessary 
for  further  research  on  this  project.  The  methods  used  to 
analyze  the  Images  from  the  AFP-675  project  should  be  similar 
to  the  methods  required  to  analyze  the  GIMI  output;  therefore, 
the  methods  covered  herein  should  support  GIMI  as  well  as  AFP- 
675.  (Carruthers,  1992) 

C.  THESIS  OUTLIME  AMD  OBJECTIVES 

In  addition  to  the  purpose  stated  above,  this  paper  is 
intended  to  provide  necessary  background  information  for 
follow  on  work.  The  paper  is  organized  as  follows. 

Chapter  I :  This  chapter  provides  general  background 
information  on  the  Far  Ultraviolet  Cameras  Experiment  and 
future  applications  of  work  associated  with  the  project. 

Chapter  II:  This  chapter  provides  background  information 
on  the  structure  of  the  ionosphere,  chemical  and  dynamical 
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processes  in  the  ionosphere,  and  a  brief  introduction  to  the 
night  airglov  phoneaena.  Topics  covered  include  diffusion 
aechanisB,  density  scale  height,  ionization  rates,  and 
photocheaical  equilibrium. 

Chapter  III:  Provides  a  detailed  description  of  the  Far 
UV  Cameras  and  the  methods  used  to  reduce  the  image  data  into 
data  files. 

Chapter  IV:  Covers  data  analysis  performed  as  part  of 
this  thesis.  Includes  a  description  of  the  means  by  which 
shuttle  position  and  attitude  were  determined  and  details  the 
programs  used  to  Interpret  the  image  data  files.  This  section 
includes  graphs  of  intensity  versus  observation  zenith  angle 
obtained  from  the  images. 

Chapter  V:  Contains  a  summary  of  findings  and  topics  for 
future  NFS  students  to  research. 
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ZZ.  ■aClOBODMD 

Ik.  THB  nVPIR  JkTMOSPHXSB 

1.  0«n«ral  D^aorlptlon 

111 _ T— paratura  strootura 

The  upper  atnosphere  can  be  described  in  many 
ways.  The  most  conmon  means  of  describing  the  layers  of  the 
atmosphere  is  by  the  temperature  profile.  Figure  (1)  shows  a 
theoretical  plot  of  temperature  versus  altitude  and  identifies 
the  various  layers. 


TEMPERATURE  K 


Figure  1:  Temperature  Profile  of  atmospheric 
layers . 
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In  the  thermosphere  the  sun's  radiation  in  the  vacuum 
ultraviolet  region  is  absorbed  by  atomic  and  molecular 
nitrogen  and  oxygen.  This  absorption  generates  heat,  as 
altitude  decreases  less  vacuum  ultraviolet  radiation  reaches 
the  atmosphere  resulting  in  lower  temperatures.  This  explains 
the  positive  temperature  gradient  in  the  thermosphere.  The 
temperature  in  the  thermosphere  reaches  a  minimum  at  about  85 
km  (The  mesopause) .  Below  this  point  ozone  becomes  an 
important  constituent  and  it  absorbs  longer  solar  radiation 
wavelengths  (2000  -300oA) .  From  85  km  to  ed^out  50  km  the 
increasing  density  of  ozone  causes  a  temperature  increase  with 
decreasing  altitude.  Also  absorbed  by  ozone  is  the  Lyman-a 
radiation  at  1216A  which  passes  through  the  nitrogen  and 
atomic  oxygen  to  an  altitude  of  about  65  km.  At  50  km  enough 
of  the  2000  -  300oA  and  Lyman-a  radiation  has  been  removed  to 
allow  a  positive  temperature  gradient  to  form  again  even 
though  the  density  of  ozone  continues  to  increase. 

This  thesis  will  deal  primarily  with  the 
atmosphere  above  100  km.  Therefore  our  discussions  will  be 
limited  to  the  thermosphere.  A  profile  of  thermospheric 
temperature  is  shown  in  Figure  (2).  This  profile  was 
calculated  using  the  equation: 

r-r.-  ( r.-ri2o>  (i) 
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where  T,  is  the  temperature  at  the  outer  edge  of  the  atmosphere 
and  T^2o  temperature  at  altitude  120  km,  z  is  altitude, 
and  s  is  .0255  km*\  The  figure  shows  temperature  calculated 
in  2km  increments  from  120  to  300  km.  T,  can  vary  from  lOOOK 
and  1800K  depending  on  the  solar  cycle.  1500K  is  used  in  this 
calculation.  A  value  of  200K  is  used  for  T^jo* 


tlf  1  I  —  I  I  I  ^  *  »  ■ 

iQ04estao«i«D«eiMMino 


Temperattrc  Jn  ICeVln 

Figure  2:  Temperature  vs  Altitude  from  120-300  km. 

T«  =  1500  K,  T,2o  =  200  K 

b.  Diffusion  Structure 

Another  means  of  describing  the  atmosphere  is  by 
the  mechanism  of  diffusion  which  dominates  at  a  given 
altitude.  Their  are  two  basic  mechanisms  of  diffusion  for  the 
neutral  atmosphere,  molecular  diffusion,  and  eddy  diffusion. 
Eddy  diffusion  is  the  turbulent  mixing  of  parcels  of  air  in  a 
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given  region.  This  is  the  dominant  method  in  the  troposphere, 
stratosphere,  and  mesosphere  where  the  atmosphere  is  more 
dense  and  acts  like  a  fluid.  In  the  troposphere  and 
mesosphere  strong  negative  temperature  gradients  allow  the 
easy  rise  of  warm  air  and  sinking  of  cooler  air. 
Alternatively,  the  dominant  diffusion  method  in  the 
thermosphere  is  molecular  diffusion.  In  this  region  the 
strong  positive  temperature  gradient  inhibits  the  vertical 
mixing  of  the  atmosphere.  The  one  exception  to  this  is  the 
region  between  85  km  and  120  km  where  atmospheric  tides  create 
strong  vertical  wind  shears  in  the  prevailing  horizontal 
winds.  The  stratosphere  has  a  temperature  gradient  similar  to 
that  of  the  thermosphere;  here  the  molecular  diffusion  rate 
approaches  the  eddy  diffusion  rate  but  eddy  diffusion  still 
dominates . 

As  a  result,  the  atmosphere  can  be  divided  into 
two  diffusion  regions;  the  homosphere  and  the  heterosphere. 
The  homosphere  is  the  region  below  90  km  and  has  a  nearly 
constant  mean  molecular  mass.  The  region  above  115  km  is 
called  the  heterosphere  and  is  characterized  by  the  mean 
molecular  mass  decreasing  with  altitude  due  to  varying  scale 
heights.  (Scale  height  will  be  covered  in  Chapter  III.) 
Figure  (3)  gives  a  general  indication  of  diffusion  types  and 
rates  versus  altitude,  and  illustrates  the  two  regions.  The 
figure  illustrates  the  crossover  point  at  115  km  and  shows  the 
variable  nature  of  the  eddy  diffusion  term.  Between  90  and 
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Flgurs  3:  Eddy  vs  Molecular  Diffusion  for  N2. 
(Heicklen,  pg  4) 


115  km  molecular  diffusion  and  eddy  diffusion  are  comparable 
as  illustrated  in  Figure  (3). 

o.  The  ionosphere 

The  ionosphere  is  so  named  because  of  the 
ionization  of  its  constituent  elements  by  solar  radiation. 
The  ionosphere  is  described  by  regions  of  varying  electron 
density.  The  ionoshere  is  divided  into  four  distinct  regions; 
the  D  region  which  is  lowest  in  altitude,  the  E  region,  the 
region,  and  the  F2  region  which  is  the  highest.  The  D  region 
extends  from  about  50  km  to  90  km  blending  in  with  the  lower 
region  of  the  E  region  with  no  sharp  boundary.  The  D  region 
experiences  total  solar  control.  It  builds  rapidly  in  the 
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morning,  maximizes  near  local  noon,  and  rapidly  diminishes  at 
siinset.  The  E  region  extends  from  the  top  of  the  D  region  to 
the  bottom  of  the  F  region,  approximately  90  '  130  km.  The 
electron  density  of  the  E  region  peaks  at  about  110  km.  The 
E  region  is  very  predictable  reaching  a  maximum  slightly  later 
in  the  day  than  the  D  region.  The  E  region  does  not  disappear 
at  night  but  does  reach  a  minimum  density.  The  F  region  is 
the  highest  ionospheric  region.  The  F  region  extends  upward 
from  about  130  km  and  has  no  clearly  defined  upper  limit. 
During  daylight  the  F  region  splits  into  two  distinct  layers, 
F^  and  F2.  The  F^  layer  is  the  lower  and  less  dense  layer.  It 
is  controlled  directly  by  solar  activity  and  disappears  at 
night  blending  into  a  single  F  layer.  The  F2  region  is  the 
highest  and  most  dense  layer  with  a  daytime  peak  altitude  of 
approximately  300  km.  The  F2  region  is  the  most  variable  of 
all  the  ionospheric  regions.  The  F  region  reaches  maximum 
density  a  few  hours  later  than  the  D  or  E  regions  and  does  not 
disappear  at  night.  This  is  due  to  the  substantially  lower 
loss  rate  of  the  constituent  ions  in  the  F  region  compared  to 
those  in  the  D  and  E  regions.  A  graph  of  electron  densities 
is  given  in  Figure  (4) . 
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Densities. 


2.  The  Mass  Continuity  Equation 

There  are  several  mechanisms  for  the  production  and 
loss  of  the  constituent  species  of  the  ionosphere;  diffusion, 
chemical  reactions,  and  photoionization  are  the  major 
mechanisms.  The  mass  continuity  equation  describes  production 
and  loss  involving  these  mechanisms.  The  one-dimensional 
approximation  to  the  mass  continuity  equation  is: 


dt  az 


(2) 
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where  dn/dt  is  the  rate  of  change  in  niusber  density,  P  is  the 
chemical  production  rate,  L  is  the  chemical  loss  rate,  and  n 
is  the  density  of  the  specific  species,  and  dT/dz  represents 
loss  due  to  diffusion.  Each  term  will  be  discussed  in  detail 
beginning  with  diffusion. 

a.  Diffusive  Bouilibrlgm 

The  major  neutral  species  (O,  Oj,  N2}  are  primarily 
affected  by  diffusion.  Chemical  production  and  loss  also 
occurs  with  these  species;  however,  their  rates  are  so  small 
in  comparison  to  the  rate  of  diffusion  that  they  can  be 
ignored  (P-nLcdT/dz) .  A  state  of  diffusive  equilibrium  exists 
when  diffusion  dominates  and  dn/dt  »  0.  This  results  in 
equation  (2)  reducing  to: 

0— H  '»> 

From  this  equation  it  can  be  shown: 

n(z)  -jj^,exp  " 

where  n(z)  is  the  number  density  of  the  species  and  has  units 
of  cm‘^.  Z  is  the  altitude  at  which  the  density  is  calculated 
and  Zq  is  the  lower  altitude  limit  (120  km  for  this  paper) . 
H*  is  the  density  scale  height  and  is  given  by: 
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^+(-i)  (^) 

JcT  T  dz 


(S) 


vhere  ■  is  the  atoaic  aass  of  the  atoa  or  aolecule,  g  is  the 
acceleration  due  to  gravity,  k  is  Boltzaann's  constant,  T  is 
temperature  (Kelvin) ,  and  dT/dZ  is  the  temperature  gradient. 
Scale  height  is  the  difference  in  altitude  that  it  takes  the 
density  of  a  species  to  drop  by  1/e. 


tU _ Photochemical  Equilibrium 

For  the  major  ionic  (N2^,  O2*,  O*,  VO*) ,  and  the 
minor  neutral  (NO,  N(^S),  and  N(^D))  species  in  the  ionosphere, 
chemical  reaction  rates  greatly  exceed  diffusion  rates; 
therefore,  dT/dz  «P-nL  and  diffusion  is  neglected.  A  state  of 
photochemical  equilibrium  exists  when  chemistiry  dominates  and 
dn/dt»0.  Thus  chemical  production  and  loss  rates  are  equal 
(P-nL^o)  and  equation  (2)  becomes: 

«) 


_ Photoioniiatioa 

Photoionization  by  solar  radiation  is  the  primary 
means  of  ion  production  in  the  ionosphere.  The  highest  rate 
of  ion  production  results  from  the  excitation  of  atoms  and 
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■ol«cules  by  solar  radiation.  Tbs  ionization  rate  (Q)  of  a 
given  species  is: 

C>(z.e)-ii(z)  Ji(z,e)Al  (7) 


where  n(z)  is  the  nuaber  density  given  above,  0|(1)  is  the 
ionization  cross  section,  and  is  the  intensity  of  the  solar 
radiation,  at  the  wavelength  l,at  the  altitude  z,  when  the 
solar  zenith  angle  is  6.  Each  of  these  iteas  is  described  in 
detail  in  the  following  paragraphs. 

The  intensity  can  be  calculated  using  the 
following  equation: 

Jj,(z,e)-Jji{-)exp-*<*'*>  (•) 

where  I|(<»)  is  the  intensity  at  the  top  of  the  ataosphere  and 
r(z,6)  is  the  optical  depth.  I|(a>)  is  tabulated. 

The  optical  depth  can  be  calculated  using  the 
following  equation: 

,(X)N{z,9)  (») 

where  a,(l)  is  the  species  absorption  cross  section  at  the 
wavelength  X  and  N(z,6)  is  the  slant  coluan  density.  The 
optical  depth  is  a  aeasure  of  the  attenuation  of  solar 
irradiance  by  the  Earth's  atmosphere.  Optical  depth  is 
dependent  on  the  specific^Xxavelength,  the  absorption  cross 
section  of  the  given  species,  the  altitude  dependent 
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concentration  of  the  given  apecies,  and  the  solar  xenith 
angle.  Where  optical  depth  reaches  unity  (t«1) ,  a  noticeable 
drop  in  intensity  (Ij(Z))  occurs.  This  drop  is  ill\istrated  in 
Figure (5) . 


Figure  5:  optical  Depth  vs  Intensity. 


The  column  density,  N(z,e),  of  the  species  and  can  be 
calculated  from  the  number  density  using  the  following 
equation: 


N(z,d)~f  n(z)_^£g. 
Jz  coed 


(10) 
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wher*  6  is  the  angle  between  the  direction  of  the  solar 
radiation  and  the  vertical  axis,  generally  called  the  solar 
zenith  angle. 

3.  Cbealoal  Produotloa  sad  Loss  ia  the  Zoaesphere 

Using  the  ionization  rate  described  above,  the  mass 
continuity  equation,  and  using  tabulated  equations  for  the 
production  and  loss  of  ions  in  the  ionosphere,  it  is  possible 
to  calculate  the  density  of  each  ionic  species.  The  constants 
associated  with  these  equations  for  production  and  loss  are 
given  in  Appendix  A.  Equations  for  determining  the  density  of 
the  major  (N2^,  ,  O* ,  NO^)  and  the  minor  (  NO,  N(^S),  and 

N(^D))  constituents  are  given  below.  The  reaction 
coefficients  are  listed  in  Appendix  A. 

The  following  reactions  produce  O^: 

0+hv  -  0*+e  Cii) 

JVa+O  -  0**N2  (12) 

The  following  eq[uations  account  for  the  loss  of  O^: 

0*+N2  -  NO**N(*S)  (13) 
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0*+Oj  -  0^*0  (14) 

Using  equation  (8)  and  reactions  (11)  through  (14),  the 
density  of  can  be  calculated  using  the  following  equation: 

^ (o> 3  t O]  (15) 

where  <2  is  the  ionization  rate  discussed  above.  The  subscript 
following  the  reaction  coefficient  K  is  the  reaction  nunber 
listed  in  Appendix  A,  and  [X]  is  the  number  density  of  species 
X.  Similar  calculations  can  be  done  for  NO*,  NO,  N(^S),  N(^D), 
Nj*,  and  Oj*. 

NO*  is  produced  by  the  following  reactions: 


N2+O  ■*  NO**Ni^D) 

(14) 

0*+N2  -  NO**N{*S) 

(17) 

p2*+NO  -  NO**N{*S) 

(18) 

02*N{*S)  -  N0*+0 

(19) 

The  following  reactions  lead  to  the  loss  of  NO*: 

NO*+e  -  0+m*S)  (20) 
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N0**e  -  o*m^D) 


(21) 


Applying  raactions  (16)  through  (21),  the  density  of  NO^  can 
be  calculated  froe  the  following  equation: 


[iWD*] 


[01  (0*1  [SOI  ♦iq,  [Q»*]  [W(*S3 


(22) 


The  following  reactions  account  for  NO  production: 


N(*S)  +O2  -  NO+O 

(23) 

NOD)  +02  -*  NO*0 

(24) 

The  reactions  for  NO  loss  are: 


o^*NO  -  m**02 

(25) 

N(*S)  *N0  -  N2*0 

(26) 

N(^D)  *N0  ■*  N2+O 

(27) 

Combining  reactions  (23)  through  (27)  yields  the  following 
calculation  for  the  number  density  of  NO: 


[M?]  -  ^8  ]  W  1  W 

K22  IC^I  +J^3  [N{*S)  ]  [N(^D) 


(28) 


18 


I',-"*, 


Thm  r«actions  for  M(*S)  production  aro: 


0  ♦Na  -  NO**N(*S) 

(29) 

liO+e  -  0*N(*S) 

(30) 

N{^D)  *e  -  m*S)  *e 

(31) 

m^D)  *o  ■*  m*s)  +0 

(32) 

N(^S)  loss  results  from  the  following  reactions: 

O  +N(*S)  -  N0*+0  (33) 


N(*S)  +O2  -  MHO 

(34) 

N(*S)  *NO  -  N2*0 

(35) 

Combining  these  equations  in  the  form  of  equation  (8)  yields; 

[N{*S)]  -  ^0  te]  W^D)  3  [e]  [jy(»g)  ]  [O]  ^ 

Ki,(Q;i  ♦Ki,[a,VJC„  [WOl 


N(^D)  is  produced  by  the  following  reactions: 


N2*0  ■*  NO*+N(^D) 


N:*e  -  N(^D)  *nOd) 


(37) 


(38) 
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NO**e  ^  o*m^D) 


(39) 


The  following  reactions  result  in  N(^D)loss: 


N(^D)  *02  -  N0*0 

(40) 

Ni^D)  *e  -  N{*S)  *e 

(41) 

N(^D)  *0  -  N(*S)  *0 

(42) 

Ni^D)  *NO  -  N2+0 

(43) 

Combining  equations  (37)  through  (43)  and  solving  for  [N(^D)] 
yields: 

fw(2nn  -  [el  [J^O-]  [e]  . 

+^30  [e]  +^31  [O]  [NO] 


Nj*  is  produced  by  the  following  reaction: 

N2*hv  -  N2*e 

The  following  reactions  result  in  N2*  loss: 

N2*0  -  NO*+N(^D) 

N2*0  -  0*+N2 


(45) 


(46) 

(47) 
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N2*0^  - 


(4t) 


-  N[^D)  *N{^D)  (4») 

Combining  these  reactions  and  solving  for  [N2^]  yields  the 
following  solution: 


_ i?(JO _ 

+jq8  IO3)  ♦jq,  [el 


(50) 


The  following  reactions  lead  to  the  production  of  O^i 


Oj+iiv  -  p^+e 

(51) 

W2+O2  -  O2+W2 

(52) 

0*+02  -  O2+O 

(53) 

lead  to  the  loss  of  02^: 

02*+N0  -  NO^+Oj 

(54) 

2+N(*S)  -  NO**0 

(55) 

Oj+e  -  0+0 

(56) 

Combining  reaction  (51)  through  (56)  results  in  the  following 
solution: 
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lCi2  +^3  itf(*S)  1  *K2t  [e] 

The  above  equations  are  used  in  photochesical  models 
such  as  the  ones  presented  in  Cleary  (1986) . 

B.  THB  VXOHT  AIBOLOW  VEBMOMBIIX 

The  night  airglow  spectrum  has  been  extensively  studied. 
Far  ultraviolet  measurements  of  the  night  airglow  spectra  have 
been  assembled  from  several  sources  including  Apollo  16, 
numerous  rocketbome  experiments,  and  several  space  shuttle 
and  Skylab  based  experiments.  The  nature  of  the  spectnun  is 
well  established.  For  the  purposes  of  this  paper,  the  far 
ultraviolet  night  airglow  spectrum  can  be  divided  into  two 
regions;  1000  -  160oA  and  1600  •  200oA.  This  serves  two 
purposes;  first  it  separates  the  spectrum  into  the  respective 
regions  of  coverage  for  the  two  cameras,  and  second,  it 
separates  the  region  of  primarily  oxygen  related  emissions 
from  the  region  of  primarily  NO  related  emissions. 

At  night,  the  free  ionospheric  electrons  and  ions 
recombine  resulting  in  a  reduction  of  ion  density.  During 
this  recombination  process  the  atom  or  molecule  can  emit 
photons  at  various  wavelengths.  Figure  (6)  illustrates  the 
night  airglow  spectra  between  750  and  140oA.  This  region  is 
characterized  by  the  spectra  of  O,  O*,  He,  and  H.  Of  note  are 
the  atomic  oxygen  recombination  emission  lines  at  1304A  and 
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1356A.  Th«se  ar*  th«  only  significant  oiissions  in  tlM 
filtarad  Casara  Ona  passband. 


Figura  6:  Night  Airglow  Spectra  between 
750  -  1400A  obtained  from  STP78-1  at  an 
Altitude  of  600  km.  (Rees  pp  140) 


Figure  (7)  shows  the  energy  level  diagram  for  01.  The 
transitions  which  are  detectable  by  filtered  Camera  One  cover 
the  range  from  10  to ’7.5  eV. 

The  nightglow  spectra  of  the  1600-2000A  region  is 
characterized  by  the  emissions  of  the  NO-f  band.  Figure  (8) 
illustrates  the  nightglow  spectziam  from  1600  -  295oA.  Note 
the  presence  of  two  NO  6-band  amissions  below  2000A.  These 
are  the  only  emissions  detectable  by  filtered  Camera  Two,  the 
upper  limit  of  the  camera  prevents  detection  of  NO  y-band 
emissions  and  allows  the  study  of  f-band  emissions  only. 
Figure  (9)  shows  the  energy  level  diagram  for  NO.  The 
passband  of  filtered  Camera  Two  equates  to  potentials  between 


(Rees,  Appendix  3) 

7.75  to  6.2eV.  In  Figure  (9),  the  N  and  O  atoas  approach 
along  the  a^IIcurve  to  the  state.  The  state  decays  to 
the  ground  state  emitting  the  6  bands.  The  reaction  is: 

N(*S)  +0(^P)  ♦hv  {b-bands)  (*•) 

The  alnina  for  the  C  and  X  states  do  not  have  the  saae 
intemuclear  separation.  This  leads  to  the  enlssion  of 
several  bands  in  the  transition.  (Rees,  pp. 140-146) 
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Heiltarrg  U, 


rigur«  •:  NightglowSpectrafrOTilS^ 
2950A.  Obtained  froa  the  S3>4 
Satellite.  (Rees  ppl41) 


ZZZ.  va  WMM  VW  CMfllM  UnSZMBHT 

The  Naval  Research  Laboratory's  Par  UV  Cameras  Experieent 
was  flown  aboard  STS*'39  bettreen  J4>ril  28  and  May  6  1991  as 
part  of  the  Air  Force  Prograe  675  payload.  The  purpose  of 
the  Far  UV  cameras  was  to  obtain  images  of  the  upper 
atmosphere,  the  shuttle  environment,  and  celestial  targets 
including  stars  and  diffuse  backgrounds.  The  wavelength  band 
of  the  cameras  was  lOSoA  to  200oA.  Two  cameras  were  used  to 
allow  observations  in  two  overlapping  passbands;  1050  >  160oA, 
and  1230  -  200oA.  Filters  were  used  on  some  images  to  allow 
imaging  in  two  distinct  passbands;  1230ol600A,  and  1650~2000A. 

A.  TEE  PAR  UV  CAMERAS 

The  electrographic  Schmidt  cameras  developed  at  NRL  are 
designed  to  provide  for  maximum  sensitivity  for  detecting 
diffuse  sources  in  the  far  UV  range.  This  is  accomplished 
through  the  use  of  (a)  low  f-number  optics,  (b)  high 
efficiency  opaque  alkali  halide  photocathodes,  and  (c) 
microchannel  intensification  to  permit  detection  of  single 
photoelectron  events.  The  objective  was  to  develop  from 
these  far  UV  images  a  method  by  which  an  altitude  profile  of 
the  upper  atmosphere  can  be  determined.  Both  cameras  have 
aperture  diameters  of  75mm  and  a  circular  field-of-view  of  20" 
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in  diaMtar.  Tha  opaqua  alkali  halida  photocathodas  hava  auch 
highar  quantua  afficiancias  in  tha  vacuua  ultraviolat  than  do 
coaaon  saaitranaparant  photocathodas.  As  a  rasult,  tha  two 
caaaras  hava  axtraaaly  high  ovarall  dataction  afficiancias. 
Tha  diffaranca  batwaan  tha  two  caaaras  lias  in  thair  diffusa 
sourca  dataction  capability  and  thair  oparating  wavalangth 
ragion. 

1.  Caaara  Oaa 

Caaara  one  covers  the  1050  -  1600  Angstroa  range. 
This  camera  uses  a  LiF  Schmidt  corrector  plate  and  a  KBr 
photocathode.  The  wavelength  range  can  be  limited  to  1230- 
160oA  by  using  a  CaFj  filter.  This  permits  the  removal  of  the 
strong  Lyman-a  emissions  at  1216A  from  the  image  thereby 
preventing  saturation  of  the  image  by  this  emission.  Camera 
One  do>3s  not  have  a  microchannel  plate  (MCP)  intensifier  and 
therefore  single  photoelectron  events  are  not  detected  by  this 
camera.  An  MCP  is  not  useful  at  this  wavelength  region 
because  the  intensity  of  Lyman-a  nightglow  is  on  the  order  of 
2000  Rayleighs  and  an  MCP  would  saturate  unfiltered  images. 
The  diffuse  source  sensitivity  of  ceunera  one  expressed  as 
units  of  optical  density  (on  the  processed  film)  per 
KiloRayleigh-second  is  estimated  at  0.04  d/kR-sec  at  a 
wavelength  of  1216A. 


.  '  . . . --'I  -  N- L,.i 

2.  CaMra  Two 

Th«  Mcond  caMra  covara  tha  1230-2000A  ranga.  This 
caaara  usas  a  LiF2  Schaldt  corractor  and  a  Csl  photocathoda. 
Tha  ranga  can  ba  lisltad  to  1650-2 OOoA  using  a  Si02  filtar. 
Photocathoda  dataction  afficiancy  varsus  wavalangth  is  givan 
in  Figura  (10).  An  (MCP)  is  usad  in  Casara  Two.  This  is 
raquirad  in  ordar  to  intansify  tha  night  airglov  of 
wavalangths  longar  than  123oA  whara  tha  typical  intensity  is 
less  than  10  Rayleighs.  The  diffuse  source  sensitivity  for 
cemera  two  is  estimated  at  4.0  d/kR-sec  at  1304A  Camera  Two 
can  detect  a  IR  emission  of  1304A  light  at  an  average  density 
of  .1  using  an  exposure  time  of  25  seconds.  The  longest 
exposure  time  was  100  seconds,  therefore  an  emission  intensity 
of  .25R  would  be  detectable  by  camera  two.  The  minimum 
intensity  of  the  far  UV  diffuse  galactic  background  is  ed>out 
2.5R  when  integrated  over  the  camera  two  passband.  The 
minimum  intensity  of  the  night  recombination  airglow 
(1356A)  is  about  IR.  As  a  result  of  the  MCP  intensification, 
single  photoelectron  events  are  observable  by  camera  two. 
(Carruthers, 1992) 
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Figur*  10:  Photocathode  Detaction  Efficiency  vs 
Wavelength. 

B.  DATA  RBOUCTXOK 

Nearly  400  images  were  recorded  by  the  two  cameras  during 
the  STS-39  mission.  After  flight  the  images  were  scanned  into 
digital  format  using  a  microdensitometer.  This  instrument 
digitises  values  of  optical  density  <i-log(J^J)  idtere  and 

I  are  the  incident  and  transmitted  intensities  of  the  scanning 
beam.  The  images  were  stored  at  NPS  on  a  Sun  Microsystem 
computer  in  the  Flexible  Image  Transport  System  (FITS)  format. 
The  tapes  are  manipulated  using  Interactive  Data  Language 
(IDL) . 
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XT.  lATlk  AK&LTSZS 

Jk.  WOMMMT  MMD  lOTBOM 

1.  Sliuttl*  Position  Dots  Piloa 

The  shuttle  position  data  is  contained  on  the 
Post  flight  Attitude  Trajectory  History  (PATH)  data  tapes.  The 
PATH  data  tapes  used  for  this  project  are  held  at  Naval 
Research  Lab  Washington  DC,  Space  Science  Division.  The 
contents  of  the  PATH  data  tapes  are  described  in  detail  in 
NASA  publication  ICD-I-TOP-002  dated  15  January  1988.  The 
required  material  from  the  tapes  was  obtained  in  November 
1992.  All  required  PATH  data  for  each  exposure  taken  during 
the  PF03  block  of  the  FAR  UV  cameras  experiment  is  contained 
in  each  record.  A  PATH  data  record  is  recorded  for  every  10 
seconds  of  flight  except  when  transmission  problems  prevented 
reception.  Although  some  breaks  in  data  were  experienced,  the 
PATH  tapes  provide  enough  information  to  determine  shuttle 
position  and  attitude  accurately  for  each  airglow  image.  A 
description  of  the  data  fields  contained  in  each  PATH  record 
is  contained  in  Appendix  B. 

Aj _ Positional  Analysis 

The  determination  of  where  the  center  of  each 
image  is  pointed  is  critical  to  this  project.  Specifically, 
the  angle  between  the  local  vertical  and  the  pointing 
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direction  of  th«  caMra  must  be  known.  Froa  this  the 
observation  senith  angle  (OZA)  can  be  detemined  for  any  point 
in  the  field-of-view.  Therefore,  the  position  and  attitude  of 
the  orbiter  aust  be  deterained.  Two  separate  coordinate 
systeas  are  used  to  define  the  position  and  attitude  of  the 
shuttle.  The  Ariesotrueoof-date  systea  shown  in  Figure (11), 
and  the  orbiter  body  axis  coordinate  systea.  The  Aries^true- 
of-date  system  is  based  on  aligning  the  -»-x  axis  with  the  right 
ascension  of  the  true-of-date  equinox  of  Aries.  The  xy-plane 
is  aligned  with  the  true-of-date  celestial  equator.  The  +z 
axis  is  upward  through  true  north,  and  the  y-axis  completes  a 
right  hand  system.  Right  ascension  is  measured  0-360  degrees 
from  the  +x  axis  positive  in  the  direction  of  the  -fy  axis. 
Declination  is  the  angle  between  the  xy-plane  and  the  radial 
vector  to  the  orbiter  and  measured  positive  toward  the  +z  axis 
and  negative  toward  the  -z  axis.  The  position  of  Aries  used 
in  this  system  has  been  updated  since  the  1958  standard  was 
established  (hence,  true-of-date) .  The  SHPOS  data  field  from 
the  PATH  data  record  gives  the  orbiter  position  in  the  Aries 
true-of-date  coordinate  system.  Right  ascension  (RA)  and 
declination  (DEC)  can  be  calculated  from  the  xyz  coordinates 
given  in  SHPOS. 

The  second  system  is  the  orbiter  body-axis  system. 
This  system  is  illustrated  in  Figure  (12) .  In  this  system  the 
orbiter  x-axis  is  aligned  longitudinal  to  the  body  and  is 
positive  out  the  front  of  the  orbiter.  The  z  axis  passes 
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through  th«  center  of  gravity  and  is  directed  out  the  top  and 
bottoai  of  the  orbiter  normal  to  the  x-axis.  The  z-axis  is 
positive  out  the  bottos  of  the  orbiter.  The  shuttle  y-axis  is 
orthogonal  to  the  xz-plane  and  completes  a  right  hand  system. 
The  -t-y  axis  is  directed  outward  toward  the  right  wing. 


Figure  IX:  Aries  True  of  Date  Coordinate  System. 

The  SHXRA  and  SHXDEC  data  fields  give  the  aimpoint  of  the 
shuttle  -t-x  axis  in  terms  of  right  ascension  and  declination. 
The  SHZRA  and  SHDEC  data  fields  give  the  aimpoint  of  the 
shuttle  -z  axis  in  terms  of  right  ascension  and  declination. 

Since  the  field-of-view  is  aligned  with  the 
orbiter  -z  axis,  the  orientation  of  the  orbiter  -z  axis  must 
be  determined.  The  first  step  is  to  determine  the  orientation 
of  the  z-axis  with  respect  to  the  orbiter  radial  (r) .  Using 
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Coordinate  Systen. 

the  unit  vector  defined  by  the  RA  and  DEC  of  the  orblter,  and 
the  unit  vector  defined  by  the  RA  and  DEC  of  the  shuttle  -z 
axis;  the  dot  product  of  the  two  conflms  the  two  vectors  are 
orthogonal.  This  proves  that  the  z-axls  Is  orthogonal  to  the 
outward  radial  from  the  center  of  the  earth.  Using  the 
Shuttle  axis  RA  and  DEC,  the  dot  product  of  this  vector  and 
the  radial  confirm  the  two  are  parallel.  Since  both  points 
pass  through  the  shuttle  position  and  are  parallel  the  shuttle 
'•■X  axis  passes  through  the  center  of  the  earth.  This  verifies 
that  the  orblter  Is  nose  down  as  expected.  The  final  question 
Is  the  position  of  the  orblter  ~z  axis  with  respect  to  the 
Arles-true»of-date  z-axls.  To  determine  this,  the  RA  and  Dec 
of  the  orblter,  the  RA  and  DEC  of  the  orblter  -z  axis,  and  the 
origin  are  used  to  define  a  plane  containing  the  three  points. 
The  angle  between  this  plane  and  the  Arles**true*of~date  -fz 
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axis  is  the  angle  of  the  orbiter  >z  axis.  Using  a  cross 
product  of  the  two  unit  vectors  to  define  a  noraal  vector  to 
the  plane  and  the  angle  between  the  Aries-true-of-date  -t-Z  axis 
and  the  normal  was  calculated.  This  angle  +/~  90  degrees 
gives  the  angle  of  the  orbiter  z-axis  with  respect  to  the 
Aries-true-of-date  z-axis.  Calculations  confirmed  the  orbiter 
z-axis  to  be  in  the  orbital  plane  which  is  at  an  inclination 
of  57  degrees.  This  confirms  that  the  base  of  the  shuttle  was 
in  the  ram  direction  as  expected.  An  illustration  of  orbiter 
attitude  is  given  in  Figure (13) .  Having  confirmed  the 
orientation  of  the  f ield-of-view,  the  next  step  is  to  begin 
image  processing. 
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rlgur*  X3:  Orbiter  Attitude  for  Airglov  laages. 
(Carruthers,  1992) 
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2 .  Zaag*  ProoMsing 


Image  processing  is  performed  using  Interactive  Data 
Language  (IDL) .  The  images  were  received  from  Naval  Research 
Led)  in  the  FITS  format  described  above  and  were  read  using  the 
IDL  READFITS  function.  Four  programs  were  developed.  Three 
are  required  to  recover  data  from  an  image.  The  fourth 
program  is  used  to  print  output.  The  programs  are  contained 
in  Appendix  C. 

The  first  program,  UVFILBS.PRO,  is  activated  by  typing 
UVFILE8, IMAGE  at  the  IDL  prompt.  DVFILB8.PR0  is  the  procedure- 
name  and  the  parameter  IMAGE  is  the  varieUsle  name  of  the 
image.  This  parameter  may  have  any  name  as  long  as  it  is 
consistent  for  all  four  procedures.  This  procedure  allows  the 
user  to  access  the  image  using  the  IDL  PICKFILE  widget.  Once 
selected,  the  image  is  read  using  the  READFITS  function  and  is 
rebinned  from  2048  x  2048  to  512  x  512  to  facilitate 
processing.  The  image  array  is  divided  by  800  to  convert  from 
photodensitometer  output  to  density  units.  A  512  x  512  window 
is  opened  and  the  image  is  presented  to  the  user  in  black  and 
white.  The  screen  will  prompt  the  user  to  select  a  color 
table.  Color  table  (13)  (Rainbow)  is  recommended.  For 
exposures  of  30  seconds,  color  table  (1)  will  also  give  a  good 
visual  presentation  of  the  image. 

Once  the  first  program  has  run,  the  second  program  is 
called  by  typing  DVBACKGROUMD, IMAGE.  This  program  allows  the 
user  to  remove  background  noise  form  the  image.  Noise  can  be 
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imparted  to  the  image  file  by  the  PDS  procedure  described  in 
the  data  reduction  section,  by  copying  the  files  over  a  data 
net,  and  background  effects  from  the  film  that  was  used.  The 
UVBACKGROUND  procedure  allows  the  user  to  define  a  box  by 
selecting  two  points  with  the  mouse  cursor.  This  is  done  in 
one  of  the  four  comers  outside  of  the  circular  image.  Care 
must  be  taken  not  to  include  any  portion  of  the  circular  image 
in  the  box  and  not  to  include  any  artifact  or  defect  which 
will  skew  results.  The  mean  density  unit  value  of  this  region 
is  calculatea  and  this  value  is  siibtracted  from  the  image 
array.  The  new  image  replaces  the  old  image  in  the  view 
window. 

The  remaining  two  progreuns  are  called  by  typing  either 
nvcURSOR, IMAGE  or  UVCUR80R1 , IMAGE .  These  programs  produce 
plots  of  intensity  verses  observation  zenith  angle.  The  only 
difference  between  the  two  programs  is  that  UVCURSORl  produces 
postscript  files  of  the  output  graphs  for  printing.  Both 
procedures  will  produce  on  screen  graphs.  After  calling  the 
program  the  user  uses  the  mouse  cursor  left  button  to  select 
points  on  the  image.  The  points  will  be  connected  by  a  line 
to  aid  the  user.  Up  to  512  points  may  be  selected.  When 
completed,  press  the  right  button  to  exit  the  cursor  loop. 
The  program  then  performs  calibration  calculations.  The 
program  asks  for  the  number  of  seconds  of  the  exposure  and  for 
calibration  factors  in  Rayleighs  and  Photons.  Calibration 
factors  are  listed  in  Table  (1).  The  program  then  calculates 
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calibrated  intensity  values  by  Multiplying  the  density  units 
for  each  point  selected  by  the  cursor  by  the  calibration 
factor  and  dividing  by  exposure  tise  in  seconds  The  UVCURSOR 
procedures  produce  plots  of  intensity  versus  observation 
zenith  angle  (OZA) .  OZA  is  based  on  the  aiiqpoint  of  the 
casera  having  an  OZA  of  90  degrees.  The  algoritha  detemlnes 
the  number  of  degrees  per  pixel  based  on  a  20  degree  field-of- 
view  and  subtracts  from  a  lower  altitude  limit  of  100  degrees 
OZA.  For  future  research,  a  correction  factor  to  account  for 
the  caunera  almpoint  being  not  equal  to  90  degrees  can  easily 
be  added  to  the  program.  For  the  benefit  of  the  user.  Two 
plots  are  produced;  the  first  gives  intensity  in  Rayleighs, 
the  second  plot  is  in  Photons/sec/steradian/cm^'  The  UVCURSOR 
procedure  can  be  run 

Table  I:  FAR  UV  CAMERA  CALIBRATION 

FACTORS.  R  -  RAYLEIGHS,  PH  -  PHOTONS. 


Camera 

Filter 

Factor 

1 

no 

3.596168E4  R 

1 

no 

2.858751E9  Ph 

1 

yes 

6.263172E4  R 

1 

yes 

4.978868E9  Ph 

2 

no 

7.074304E2  R 

2 

no 

6.623675E7  Ph 

2 

yes 

4.444141E3  R 

2 

yes 

3.532842E8  Ph 
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multiple  times  hovever  the  two  plot  windows  must  be  quit 
first.  After  running  the  UVCURSORl  procedure  the  graphics 
driver  is  set  to  postscript,  this  must  be  reset  by  the  user  to 
the  default  driver  for  his  system  before  additional  image  work 
can  be  done.  Therefore,  the  UVCURSORl  procedure  should  be  the 
last  one  used  for  a  given  image,  and  only  once  the  user  is 
ready  for  hardcopy  plots.  One  additional  plot  is  available; 
a  contour  plot.  This  plot  is  produced  by  the  procedure 
UVCOMTOUR  and  is  called  by  typing  UVOOMTOUR^ZIIAOE.  This  plot 
is  in  Rayleighs.  As  with  UVCURSORl,  this  procedure  should  be 
run  after  UVBACXGROUllD  and  the  graphics  driver  must  be  reset 
after  the  procedure  is  run. 

B.  RESULTS 

Six  airglow  images  were  obtained  from  NRL,  all  are  filterd 
Camera  Two  images.  The  files  are  summarized  in  Table  (II). 

Table  II:  FAR  UV  FILENAMES,  EXPOSURE  LENGTHS,  AND  TIME  OF 

IMAGING. 


Filename 

Exposure  Length 

GMT 

FARUV462 . IFITS 

3  seconds 

05.04.16.21.59 

FARUV143 . IFITS 

10  seconds 

05.04.16.22.09 

FARUV144 . IFITS 

30  seconds 

05.04.16.22.39 

FARUV62 . IFITS 

3  seconds 

05.04.16.23.11 

FARUV63 . IFITS 

10  seconds 

05.04.16.23.21 

FARUV55. IFITS 

30  seconds 

05.04.16.23.51 
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K  color  display  of  thssas  iBagss  is  produced  using 
UVrZLBSfZiaMII.  Figure  (14)  shows  a  contour  plot  of  onr  of 
these  images  (FARUV144.IFITS) .  This  plot  was  generated  using 
UFOOMTOUR, ZllhGB .  Next,  the  orientation  of  the  orbiter's  *i'X 
and  ~z  axes  were  established.  From  this  attitude  information 
a  method  for  determining  intensity  versus  observation  zenith 
angle  was  developed.  A  typical  exaiq;>le  of  this  is  shown  in 
Figure  (15) .  This  figure  shows  the  intensity  profile 
corresponding  to  the  rectangular  area  in  Figure  ( 14 ) .  This 
profile  and  printout  was  produced  using  UVCTOSOR/ZMAOB. 
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Contour  Plot  of  For  UV  Airqlow  Intensity 
(Level  X  100) 


0  1 00  200  300  400  500  600 


Figure  14:  Contour  Plot  of  Intensity  vs  Bin.  Rectangular 
Area  shows  the  Vertical  Coluaun  froa  which  the  Plot  in 
Figure  (15)  was  derived. 
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Observation  Zenith  Angle  in  Degrees 
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Figur*  15;  Intensity  vs  OZA  for  the  Area  shown  in  Figure 
(14)  above. 
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A.  SiniiaUlY  OF  FZXDZMOf 

Shuttl*  position  and  oriantation  has  bsan  sstablishsd  for 
six  isagas  of  tha  Far  Ultraviolat  Caaaras  Expariaant.  All 
iaagas  ara  filtarad  Caaara  Tvo  iaagas.  Tha  oriantation  of  tha 
orbitar's  -t-x  and  •z  axas  wara  astablishad  and  froa  this  a  plot 
of  intensity  versus  observation  sanith  angle  was  produced. 

Froa  these  intensity  plots,  accurate  aaasureaents  of 
thermospheric  and  ionospheric  densities  can  be  derived.  From 
these  densities  an  approximation  of  free  electron  density  can 
be  made.  Since  there  are  negligible  numbers  of  negative  ions; 
charge  equilibrium  requires  that  the  number  of  positive  ions 
equal  the  number  of  free  electrons.  Therefore,  tha  density  of 
positive  ions  equals  the  density  of  electrons.  From  this,  a 
prediction  of  RF  frequency  effects  can  be  made. 

B.  TOPICS  FOR  FURTHER  STUDY 

Since  this  is  the  first  NPS  thesis  on  the  Far  Ultraviolet 
Cameras  Experiment,  there  are  numerous  topics  for  further 
research. 

1.  Determination  of  thermospheric/ ionospheric  densities 
from  the  intensity  plots  and  relating  them  to  geographic 
regions.  This  would  include  the  removal  of  background 
radiation  from  celestial  sources. 
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2.  Deriving  an  algoritha  for  procaasing  tha  data  fr<»  tha 
Global  laaging  Monitor  of  tha  lonosi^ara  aatallita  basad  upon 
tha  Par  UV  Caaera  inagas.  Usa  of  this  algor itha  to  provide 
real  tiae  RF  predictions  to  a  geographic  region. 

3.  Spatial  variations  in  charged  particle  densities  in 
the  ionosphere  can  be  detemined  froa  spatial  inforaation  in 
the  iaages. 

4.  Any  topics  concerning  the  aurora  using  six  additional 
iaages  of  the  aurora  obtained  fora  NRL.  This  will  require 
additional  position  analysis  as  all  aurora  iaages  were 
giabled . 

5.  Additional  unfiltered  iaages  are  availedale  for  work 
concerning  Lyman-a  effects  and  01. 
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Chaaical  Raaction  Rates 


This  section  contains  all  the  cheaical  reactions  and  reaction 
rate  coefficients  used  for  the  density  calculation  prograas 
for  the  neutral  ataosphere  aodel.  Sources  for  the  equations 
are  given  on  the  last  page  of  the  appendix.  Equation  nuabers 
correspond  to  those  in  Cleary  (1986) . 
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Ch«Bical  Reaction  Rates  (cont) 
Nueber  Reaction  Rate  (oa’  s*^) 


1 

0  +  hv  -»  0^+  e 

I,  «  9.1E-7  s  ’ 

2 

I,  «  1.6E>6  s  ’ 

3 

Nj  +  hv  -  +  e 

I,  -  1.2E-6  s*’ 

17a 

N2*  +  0  -►  NO*  +  N(2d) 

1.4E-10  (T/300)*®  ** 

17b 

Nj*  +  0  -»  0*  +  Nj 

l.OE-11  (T/300)®“ 

18 

5.1E-11  (T/300)®® 

19 

Nj*  +  e  -  N(2d)  +  N(2d) 

1.8E-7  (T/300)®” 

20 

0*  +  N2  -*  NO*  +  N(*S) 

5.0E-13 

21 

0*  +  O2  -*  O2*  +  0 

2.0E-11  (T/300)'®* 

22 

4.4E-10 
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Cheaiical  Reaction  Rates  (cont) 


NuBber 

Reaction 

Rate  (CB^  8*^) 

23 

0^*  +  N(*S)  NO^  +  0 

1.8E-10 

24 

O2  +  e  -*  0  +  0 

l.OE-5 

26 

NO*  +  e  0  +  N(*S) 

.25M.2E-7  (T/300)'®" 

27 

NO*  +  e  -»  0  +  N(*D) 

.75*4.2E-7  (T/300)*®-*® 

28 

N(*S)  +  Oj  -»  NO  +  0 

4.4E-12exp(-3220/T) 

29 

N(*D)  +  Oj  NO  +  0 

6.0e-12 

30 

N(2d)  +  e  -»  N(*S)  +  e 

6.0e-10  (T/300)®-* 

31 

N(2d)  +  0  -  N(*S)  +  0 

5.0E-13 

33 

N(*S)  +  NO  -  Nj  +  0 

3.4E-11 

34 

N(^D)  +  NO  -►  Nj  +  0 

7.0E-11 

Sources  of  Cheaical  Reaction  Rates 


Equation  nunber  and  source 

1.  Torr  et  al  (1979) 

2.  Torr  et  al  (1979) 

3.  Torr  et  al  (1979) 

17a.  McFarland  et  al  (1974) 

17b.  McFarland  et  al  (1974) 

18.  Lindinger  et  al  (1974) 

19.  Mehr  &  Biondi  (1969) 

20.  Lindinger  et  al  (1974) 

21.  McFarland  et  al  (1973) 

22.  McFarland  et  al  (1974) 

23.  Goldan  et  al  (1966) 

24.  Torr  et  al  (1976) 

26.  Torr  et  al  (1976)  Kley  et  al 

27.  Torr  et  al  (1976)  Kley  et  al 

28.  Becker  et  al  (1969) 

29.  Lin  &  Kaufman  (1971) 

30.  Frederick  &  Rusch  (1977) 

31.  Rusch  &  Sharp  (1981) 

33.  Lee  et  al  (1978) 

34.  Black  et  al  (1969) 


(1977) 

(1977) 
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APPnroZZ  Bt  PBTH  DBTA  FZBLD8 


This  appendix  contains  a  description  of  the  data  fields 
used  in  the  Post flight  Attitude  Trajectory  History.  A  sore 
complete  description  of  the  entire  PATH  product  can  be  found 
in  the  PATH  product  documentation  provided  by  NASA. 


Record  ntunber 

Record 

Orbit 

Orbit  number 

GMT 

YYMMDD.HHMMSS 

ATTMAT 

3x3  Transformation 
matrix. 

SHPOS 

Shuttle  position  in  the 
Aries  true  of  date 
coordinate  system. 

SHVEL 

Shuttle  Velocity 
components 

SHXRA 

Shuttle  X  axis  right 
ascension 

SHXDEC 

Shuttle  X  axis 
declination 

SHZRA 

Shuttle  +Z  axis  right 
ascension 

SHZDEC 

Shuttle  +Z  axis 
declination 

SHLON 

Shuttle  Longitude 

SHLAT 

Shuttle  Latitude 

SHZSUN 

Angle  between  shuttle  +Z 
axis  and  Sun 
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SMZEAR 

L$SNVS 

L$BBEC 

TANRAO 


Angle  between  shuttle  *Z 
axle  and  center  of  Earth 

0  >  Sun  not  visible 
1  -  S\m  visible 

0  >  Not  in  Electrical 
contasination  band 
1  -  In  contamination  band 

Tangent  radius 
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APPBIIDIX  Ct  ZMAOB  PSOCB8SIIIO  PBOQRJUI8 


This  appendix  contains  IDL  programs  used  to  read  the  Far 
UV  images  and  create  the  plots  of  intensity  versus  observation 
zenith  angle. 


pro  uvf iles , image 
f ile*pickfile (/read) 

result-readf its (file, head) 
image*rebin (result, 512, 512) 
image-image/800 . 

window, 1, xsize-512, ysize-512 
tvscl, image 

read, (' select  desired  color  table') 
read, clrtbl 
loadct, clrtbl 

return 

end 
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.,Ji^!i|||  J,u  !  4  .l4iL|it^||p|pP^P 


pio  uvbackground, image 


cursor, f , g, /device, /down 
cursor, £1 , gl, /device, /down 

j-min ( (f , f 1] ) 
jl-max{ (f , fl] ) 

k-min ( (g, gl] ) 
kl-max ( [g, gl] ) 

cornerl-image ( j : jl, k:kl> 
areal«‘n_elements  (cornerl) 
cornerl=total (cornerl) 
cor nerl“Cornerl /areal 


image"image -cornerl 
tvscl, image 

return 

end 
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,  1  ■ 


(•I  '.v  .1  ■■ 


y' 111.  i'viiicn  Vo  solt'v't  j'v'iiH,  f'jess  right  button  at  laut  point' 


oiji;ot,x,y,  /device,  /down 
count  or “0 
a»£ltarr  (512) 
b-f Itarr  (512) 


while  (!err  ne  4)  do  begin 

a  (counter) oy 
b (counter ) -image (x, y) 

cursor, xl, yl,  /device,  /down 
plots, (x,xl], (y,yl],  /device 

counter-counter-tl 
x-xl  t  yyl 

endwhile 


a-a  (0  tcounter) 
b-b(0: counter) 

len“n_elements (a) -1 

degrees-20 . /512 . 
za-lOO-  (a*degrees) 

read, ('Enter  exposure  time  in  seconds:  '), seconds 

read, ('Enter  Calibration  Value  for  Rayleighs:  '),calfactor 

read, ('Enter  Calibration  Value  for  Photons:  '), photons 

intensity-b*calf actor /seconds 
intensityp«b*photons/ seconds 

« 

window, 0, xsize-512, ysize-512 
plot,  intensity,  za,psyni-l,  $ 

xrange- (min (intensity) , max (Intensity) ] ,yran9e-[max(za) ,min (za) ] , $ 
title-' Intensity  vs  Zenith  angle', $ 
xtitle-'  Intensity  cRayleighs/sec/  steradian/cm''2 ' ,  $ 
ytitle-'Observation  Zenith  Angle  in  Degrees' 


window, 2, xsize-512, ysize-512 
plot , intensityp, za, psym-1 , $ 

xrange- (min (intensityp) , max (intensityp) ] ,  yrange-(max (za) ,min (za) ] , $ 
title-' Intensity  vs  Zenith  angle', $ 
xtitle-'  Intensity  !cPhotons/sec/steradian/cm'‘2' ,  $ 
ytitle-'Observation  Zenith  Angle  in  Degrees' 


return 

end 
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I  <•  »  t,  ( 


t'Ut  tori  » s«-l<‘ct  p^'int,  pr^ss  right  t'Ut  l  on  at  list  j 


'  ’!•  ••  l«»t  I  m.'ii-..- 


cuisoi.x.y,  /J^vic^,  /down 
oout'.t  01  •  0 

a-::tarr(512) 
b-f Itarr (512) 


while  (lorr  ne  4)  do  begin 

a ( counter ) -y 
b (counter) “image (x, y) 

cursor, xl, yl,  /device,  /down 
plots, (x,xll, [y,yll,  /device 

counter  “Counter'll 
x-xl  t  y-yl 

endwhiie 


a-a (0 : counter ) 
b-t (C : counter ) 

ler.-n_elements  (a)  -1 

degrees“20 . /512 . 
ra-lCO- (a’degrees) 

read, ('Enter  exposure  time  in  seconds:  '), seconds 

read, ('Enter  Calibration  Value  for  Rayleighs:  '),calfactor 

read, ('Enter  Calibration  Value  for  Photons:  '), photons 

intensity“b*cal factor /seconds 
intens it yp“b*photons/ seconds 

window, 0, xsire-512, ysi2e“512 
plot, intensity, za,psym“l, S 

xrange- [min (intensity) , max (intensity) ] ,  yrange“ [max (za) ,min ( za ) ] , $ 
title-' Intensity  vs  Zenith  angle', S 
xt itle-' Intensity ! cRayleighs/sec/steradian/ cm*2' , S 
ytitle-'Observation  Zenith  Angle  in  Degrees' 


window, 2, xsize-S12, ysize-512 
plot, intsnsityp, za,psym“l, S 

xrange- [min (intensityp) ,max (intensityp) ) ,  yrange- (max (za) ,min (za) } , $ 
title-' Intensity  vs  Zenith  angle', $ 
xtitle-' Intensity !cPhotons/sec/steradian/cm''2' ,  $ 
ytitle-'Observation  Zenith  Angle  in  Degrees' 
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S>'t  p;.-' 
dev  : 

plot , 

device, 

device, 
plot,  ;r. 


; , 'ps' 

: .  lonatre- '  inien]  .ps' ,  /port  rait 
lensityp,  za,  ps yiii-l,S 

X range* (min (intensityp) ,max (intensitypl  ] ,  y range* im3x(za), mi n (za) ] ,$ 
t It le*' Intensity  ws  Zenith  angle', S 
xtitle*' Intensity ! cPhotons/sec/steradian/cm*C' ,  $ 
ytitle*'Observation  Zenith  Angle  in  Degrees' 

close 


rilenarr*-'  inten2  .ps' ,  /portrait 
tensity, za.psym-l, S 

X range* (min (intensity) , max (intensity) ] ,  y range* (max (za) , man (za) ] , < 
title-' Intensity  vs  Zenith  angle', S 
xtitle-' Intensity !cIUyleighs/sec/stetadian/cm*2' , $ 
ytitle-'Observation  Zenith  Angle  in  Degrees' 


device,  close 


pro  uvcontour , image 


read, 'enter  exposure  time  in  seconds' , seconds 
read, 'enter  Rayleigh  calibration  factor', cal 

image 1- image * ca 1 / seconds 
image 1-smooth ( image 1, 15) 


set_plot, 'ps' 

device, filename-' uvcontour .ps' , /portrait 

contour, imagel, levels- (50, 100, 150, 200, 250, 300] , c_annotation- ( ' .5','!', '1.5', '2. O', 
title-' Contour  Plot  of  Far  UV  Airglow  Intensity ! c (Level  x  100)' 

device, /close 

return 

end 
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